Ordering of the lamellar phase under a shear flow 
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The dynamics of a system quenched into a state with 
lamellar order and subject to an uniform shear flow is solved in 
the large- N limit. The description is based on the Brazovskii 
free-energy and the evolution follows a convection-diffusion 
equation. Lamellae order preferentially with the normal along 
the vorticity direction. Typical lengths grow as 7 i 5/4 (with 
logarithmic corrections) in the flow direction and logarithmi- 
cally in the shear direction. Dynamical scaling holds in the 
two-dimensional case while it is violated in D — 3. 

PACS numbers: 64.75.+g; 05.70.Ln; 47.20.Hw 



Many systems in nature exhibit lamellar order. One 
example is a diblock copolymer melt where chains of 
type A and B covalently bonded end-to-end in pairs seg- 
regate at low temperatures with A-B junctions forming 
a stack of lamellae 0. In ternary mixtures a lamellar 
phase is stable with ordered sheets of surfactant separat- 
ing (say) oil and water domains Q . Lamellar order is also 
observed in Raleigh-Benard cells, where convective rolls 
form above the convective threshold j|. Further exam- 
ples include smectic liquid crystals Q , dipolar fluids with 
long-range interactions |5| and chemically reactive binary 
mixtures ||. A theoretical model for the general de- 
scription of the lamellar-disordered phase transition was 
proposed by Brazovskii Q who showed the first-order 
character of the transition induced by fluctuations. 

Lamellar phases under an applied shear flow show a 
very rich behavior which is relevant for many applica- 
tions A variety of transitions in morphology and 
orientation occur as shear rate and temperature are 
changed. Stable configurations of lamellae lying along 
the flow with the normals differently oriented have been 
observed (^Jl^l and analyzed evaluating the effects of the 
flow on the fluctuation spectra [ pTp^p^ , p^ |. 

Non-statio nary properties are far less considered 
In this Letter we study the effects of a 
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simple planar shear flow on the ordering of the lamel- 
lar phase in a system quenched from an initially high 
temperature disordered state. We present the first ana- 
lytical results on the kinetics of this system by solving 
the Brazovskii model in the limit of an infinite number 
of components of the order parameter. This is one of the 



few methods allowing an explicit solution for phase or- 
dering systems | fl9|| , A similar approach for fluids under 
shear flow has been used in @^,||,||,||,||||,||| . 

The behavior of quenched binary mixtures without im- 
posed flows is characterized by dynamical scaling: The 
structure factor obeys C(k,t) = R{t) D f[kR(t)] where 
R(t) ~ t a is the typical domain size and D is the space 
dimensionality (2^]nJ . In the case of a fluid with lamel- 
lar order, if the order parameter is not conserved as in 
the Swift-Hohcnbcrg model for Raleigh-Benard convec- 
tion p9|, r egimes exhibiting dynamical scaling have been 
found ||30[ . In models with conserved order parameter, 
when diffusion is the only segregating physical mecha- 
nism, simulations show the entanglement of the fluid into 
frozen intertwined structures [ [jTp^ ] . In this case the ef- 
fects of hydrodynamical modes are crucial for reaching 
order on large scales |33|. 

When shear is applied, our results show that lamel- 
lae grow preferentially with the perpendicular orienta- 
tion p^ ] , namely along the plane formed by the flow and 
the shear (velocity gradient) directions. Their typical 
size, obtained from the second momentum of the struc- 
ture factor, grows as 7£ 5 / 4 Vhii in the flow direction and 
as V^nt in the shear direction. Surprisingly we find that 
dynamical scaling is obeyed in two dimensions but not in 
D = 3. Our results concern the cases of conserved and 
not conserved order parameter, and apply to most of the 
systems mentioned above. 

We consider the Brazovskii free energy 

H4>} = /^{^ 2 + ^ 4 -^|V0| 2 + i(vw} (i) 

where <f>(f, t) is the order parameter field and u > 0. With 
a negative value of r and b < the system orders in one 
of the two minima of the local potential. However, when 
b > interfaces are favoured and a modulated state with 
wavevector ku — y/b/2 is stable (?]] . 

When a flow is imposed, the kinetics can be described 
by the convection-diffusion equation J34| 
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and r is a mobility coefficient, p = describes systems 
with non-conserved order parameter (NCOP) and cor- 
responds to the Swift- Hohenberg equation |^{|; the case 
p = 2 is for conserved order parameter (COP) and ap- 
plies for example to copolymer melts. For uniform shear 
flow in the x-direction v x = 7 y, 7 being the shear rate. 
Eq.(|]) neglects thermal fluctuations and possible effects 
due to differences in viscosities between the two compo- 
nents |12|. The complete description should take into 
account the coupling of Eq. (H) with the Navier-Stokes 
equation. However, the study of Eq. (|J) is a prerequi- 
site for any more general theory and usually applies to 
intermediate temporal regimes, as recognized for phase 
separation of simple binary mixtures |19|| . 

Eq. (||) can be studied analytically generalizing the 
field 4> to a vector order parameter with N components 
and taking the large- iV limit J35|. One obtains the fol- 
lowing equation in Fourier space for the structure factor 

C(M) = fa(MM-M)) 
dC(k,t) _ dC(k,t) _ 

7 X dky ~ 

- 2Tk p [r + uS(t) - bk 2 + k 4 ] C(k, t) (4) 

where we have dropped the component indices due 
to internal symmetry. S(t) has to be computed self- 
consistently through 



S(t) 



dk 



|fc|<A 



(27T) 
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C(k,t), 



(5) 



A being an ultraviolet cut-off. For a symmetric mixture 
quenched from a high temperature homogeneous phase 
an appropriate initial condition is C(k, 0) = A, where A 
is a constant. Eq. (ji|) can be integrated yielding 

C(k t) = Ae" 2r J tdT ^(r)[K:V)-bKV)+H-«s(t-T)] ^ 
where fC(r) = k + jTkxCy. Defining 

Q(t)= f dT[r + uS{t-T)} (7) 
Jo 

the analysis can be carried out asymptotically through 
the ansatz 
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Q(t) = k 4 M t- ^p-plni-lnlnt + lni;,) 



(8) 



where v p is a constant to be determined. Eq.(|^) will be 
justified a posteriori by proving the solution of the self- 
consistency problem. Specifically, the large time value 
Soo calculated through (jjj), inserted into (R), must satisfy 
(pi) asymptotically, namely 



Sac = ( k M - r )/ u 



We now briefly illustrate this for COP in D 
setting T = l = u = b = —r. 
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FIG. 1. The typical lengths as a function of the shear strain 
in the x (*), y (o) and z (•) directions. The lines are propor- 
tional to 7t 5/4 (logi) 1/2 and (logt) 1/2 . 

The symmetry of the problem suggests the use of 
the cylindrical set of variables {k x ,kj_,9), where k y = 
k± cos 9 and k z = k± sin 9. Naive power counting applied 
to Eq. (^) indicates that wave vectors in the flow direction 
scale asymptotically as k x ~ £~ 3 / 2 while in the perpen- 
dicular directions k± — /cm ~ t~ x l 2 . Then, introducing 
the new variables X — k x L\\{i), Q = (k± — kM)L±(i), 
with L||(t) = 2y / 2/3k 2 M jt 3 / 2 and L±(t) = 2\/2tfc M , and 
taking into account Eq. (||), Eq. (^|) reads as 

\8)+g(X,Q,8,t)] 



with 



C(k,t) = v 2 Af- e-^ x <' 
hi 



f = Q 2 + VSQXi 



X 2 cos 2 f 



(10) 



(11) 



and g(X,Q,9,t) = ^ 2 9j-> where gj are polyno- 

mials in X, Q, 9, \ogt. For long times g can be shown to 
provide a relevant contribution gi ~ y|rr 
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only in the region Q ~ 0, cos 9^0. Then the integration 
over Q and 9 can be performed yielding 



S(t) 



V3A 



V 2 
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dX< 



I (X 2 /8) (12) 



where Io(z) is a Bessel function. The integral of Eq. (p^), 
evaluated by the asymptotic expansion oflo(z) at large z, 
behaves as In tj y/n. Hence at large times S(t) approaches 
the constant Soo = ■ A comparison with Eq.(j|) 

OZ7T 

327r 2 fc^ 7 (4 f + l) 



(9) fixes V2 



V3A 



and verifies the ansazt. 
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Due to the presence of g, the function C(k,t) of 
Eq.(|lO|) cannot be cast in a scaling form. However, dif- 
ferent scaling behaviors are obeyed in the regions where 
the functions / or g can be respectively neglected. In 
particular, for Q ~ 0,cos# ~ 0, where g dominates, 
C(k, t) scales in the x-direction with respect to the length 
Z|| (t) = (2/5) 1 / 4 VfcI7 7 t 5 / 4 different from L\\ (t). In D = 2, 
on the other hand, where the contribution of the func- 
tion g is always negligible, one finds the scaling form 
C(k,t)=v 2 At 2 e-M x ^M. 
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FIG. 2. The structure factor is shown for the COP case on 
the planes k x — 0, k y — 0, fc 2 =0, k z — ku and for the 
NCOP case on the planes k y — 0, k z — kM- 



Ordering properties are usually inferred from the 
momenta of the structure factor. We define R x — 
[J dkC(k, t)/ j dkk 2 C(k, t)] 1 ^ 2 and similarly for the shear 
and vorticity directions. In the absence of dynamical 
scaling definitions of R x based on different momenta of 
C(k,t) may lead to different results. However in this 
case it can be proven that changing the order of the 
momentum does not change the growth exponent but 
only the logarithmic correction. ^From Eq. (|l^) we find 
Rx ~ 7t 5/4 Vhii, R y ~ Vhii and R z ~ k M x . The behav- 
ior of R x shows the relevance of the scaling with respect 
to Z|| (t) for the ordering of the system. The growth of R y 
indicates that lamellae order preferentially in the plane 
x-y. In D = 2 we find R x ~ yt 3/2 and R y - fc^/. The 
same results are obtained with non-conserved order pa- 
rameter. The behavior of R X1 R yi R z resulting from the 
numerical integration of Eq. (^) with an adaptive grid 
algorithm is shown in Fig. 1. After the initial isotropic 
evolution, the shear-induced anisotropy becomes evident 



for values of the strain yt larger than one and agreement 
with the analytical behavior is observed. 

It is also useful to illustrate the behavior of the struc- 
ture factor ( fL0| ) shown in Fig. 2. The maxima of C{k,t) 



in / 



are located at k 2 . — 5 2 — 



ky — 0, k z 



k 2 M for COP 



and at k 7 



0,k 2 y 



k 2 M for NCOP. The shape of 



C(k,t) on the planes k x = and k z — is qualitatively 
similar for COP and NCOP. At k x — 0, in particular, the 
wavevectors are not affected by the flow and the struc- 
ture factor has a circular structure with radius /cm- At 
k y = or k z — kM the behavior of C(k,t) depends on 
the conservation law. The maxima of C(k,t) developed 
with NCOP are splitted with COP into a pair of sym- 
metric narrowing peaks. This pattern is typical of the 
case with conserved dynamics and has been observed in 
other segregating systems under shear flow |37|. 



i o 



1 




1 



f 



FIG. 3. The excess viscosity as a function of the shear 
strain. The straight line has slope -2. 

Finally we turn to the study of the rheological prop- 
erties. In phase separation of binary mixtures the flow 
acts against the surface tension inducing stretching of 
domains followed by break-up processes and burst of 
small bubbles |^,^9|. An excess viscosity At] is mea- 
sured fjcjj : it reaches a maximum generally for yt > 1 
and later decays. In our case we measure the evolu- 
tion of the shear stress a xy deducible in a general way 
from (|l|) in terms of the structure factor; it is given by 
a xy (t) = J j0yok x k y (2k 2 -b)C(k,t) {24). An excess vis- 
cosity can be introduced as Ar/(t) = —a xy (t)/^ [^l). Us- 
ing Eq. jlp| ) we find Aij ~ 7~ 1 £~ 2 . The numerical results 
of Fig. 3 show that, after reaching a maximum, the de- 
cay of Ar/ agrees with the analytical behavior. The other 
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rheological indicators can be shown to behave similarly. 

In conclusion we have studied the ordering kinetics of 
a lamellar phase in shear flow by solving the dynamics 
of the convective-diffusion equation in the large-N limit. 
Regarding the debated question of the stable orienta- 
tion, our results are in agreement with the expectation 
of the stability of the perpendicular phase at high shear 
rates [fL2|. In D = 3 a logarithmic growth law is found 
along the shear direction while R x ~ 7i 5 / 4 v / In7. Inter- 
estingly, the same exponent 5/4 is found for simple bi- 
nary fluids |^6|,^J • Our results show that the scaling de- 
pends on dimensionality and is violated in D = 3. This is 
at variance with the case without shear where the same 
approach gives a scaling form independently of D [EgJ . 
Violation of scaling is quite uncommon p3| , |l9| . In this 
context this phenomenon may be related to the existence 
of two lengths Lu , lu growing with different exponents. 

A natural question is the relevance of the behavior 
of the large-N model to the physical case with N = 1. 
As without flow, symmetry arguments suggest an anal- 
ogy between the coarsening behavior of lamellar phases 
and that of vectorial models and one would argue that 
the large-N results give reliable indications also for the 
physical case 0]; in particular the large-N exponent is 
the one expected asymptotically |3(J . Simulations of the 
scalar case with shear could elucidate this point; how- 
ever strong finite size effects make the evaluation of the 
growth exponent difficult (see (3^] and refs. therein). We 
also mention that it will be important to study how the 
presence of hydrodynamics affects the picture provided 
in this Letter. 
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